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Ultra-light scalar fields, with masses of between m = 10~ 33 eV and m = 10 -22 eV, can affect the 
growth of structure in the Universe. We identify the different regimes in the evolution of ultra-light 
scalar fields, how they affect the expansion rate of the universe and how they affect the growth rate 
of cosmological perturbations. We find a number of interesting effects, discuss how they might arise 
in realistic scenarios of the early universe and comment on how they might be observed. 



I. INTRODUCTION 

Ultra-light scalar fields arise generically in high energy 
physics, most commonly as axions or other axion-like par- 
ticles (ALPs). They are the Pseudo-Goldstone bosons 
(PGBs) of spontaneously broken symmetries, which only 
acquire mass through non perturbative effects. In generic 
string theory compactifications we expect a number of 
axions pQ, similar to the well known QCD axion [2141 1 j . 
and more recently e.g. |12j . and their cosmology has 
been well studied in the past (see, for example, |13|). 
However their symmetry breaking scale is much higher, 
f a ~ 10 16 GeV, being set by the string scale and the in- 
stanton action. The axions arise from closed two cycles 
in the compact space: the symmetry breaking scale is 
only weakly dependent on the size of the cycle, so will 
be roughly constant across all the axions in a given com- 
pactification volume, whereas the mass of each axion de- 
pends exponentially on the size of the cycle and so we 
expect axions masses to evenly distribute on a logarith- 
mic mass scale all the way down to the Hubble scale 
today, H ~ 10" 33 eV pH]. 

We can characterise the Lagrangian of such a generic 
ALP using two parameters: the symmetry breaking scale, 
f a , and the overall scale of the potential, A. The axion 
field is an angular variable, 9, since the path integral 
is unchanged by the shift symmetry 9 — > 9 + 2n. The 
effective four dimensional Lagrangian can be written as: 

C=^-(d9) 2 -A 4 U(9) (1) 

where U{9) is some periodic potential. Bringing the ki- 
netic term into canonical form we define the field <p = fa9, 
with Lagrangian: 

C=\{d^ + V{ ( f ) ) (2) 

where V(4>) is again a periodic potential. Expanding the 
potential in powers of (j>/ f a , all the couplings of the field 
</> come suppressed by the large scale f a , and to quadratic 
order we find that the mass is given by: 



Hence, any axion is equivalently parameterised by its 
mass and symmetry breaking scale. 

Production of cosmological axions proceeds by the vac- 
uum realignment mechanism. When the Peccei-Quinn- 
like £7(1) symmetry is broken at the scale f a the axion 
acquires a vacuum expectation value, 9i, uncorrelated 
across different causal horizons. However provided infla- 
tion occurs after symmetry breaking, and with a reheat 
temperature T < / a , then the field is homogenised over 
our entire causal volume. This is the scenario we con- 
sider in this paper. The field 9 is a PGB and evolves 
according to the potential U{9) induced after symmetry 
breaking by instantons. Once the mass overcomes the 
Hubble drag the field begins to roll towards the mini- 
mum of the potential, in exact analogy to the minimum 
of the instanton potential restoring CV invariance in the 
Peccei-Quinn mechanism for the QCD axion. Coherent 
oscillations about this minimum lead to the production 
of the weakly coupled ALPs. [I5HT7] 

In this paper, because of the weak couplings caused 
by the high scale / a , we will choose to work simply with 
the quadratic part of the potential, and largely ignore 
the effects of anharmonicities, although we discuss them 
briefly in Section |Hl] In line with this choice we also make 
the generalisation to work with ultra-light scalar fields 
that do not respect the shift symmetry and therefore have 
no need for a periodic potential. 

Scalar fields with masses in the range 10 -33 eV < m < 
10~ 22 eV are also well motivated dark matter candidates 
and constitute what Hu has dubbed "fuzzy cold dark 
matter" , or FCDM [18] . The Compton wavelength of the 
particles associated to ultra-light scalar fields, in natural 
units, A c = 1/m, is of the size of galaxies or clusters of 
galaxies, an so the uncertainty principle prevents localisa- 
tion of the particles on any smaller scale. This naturally 
suppresses formation of structure and serves as a simple 
solution to the problem of "cuspy halos", and the large 
number of dwarf galaxies, which are not observed and are 
otherwise expected in the standard ACDM cosmological 
model. 

The large phase space density of ultra-light scalar fields 
causes them to form Bose-Einstein condensates (see [19] 
and references therein) and allows them to be treated 
as classical fields in a cosmological setting. This could 
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lead to many interesting, and potentially observable phe- 
nomena, such as formation of vortices in the conden- 
sate [101 HI], and black hole super radiance [H [H [53] 
which could provide direct tests of the "string axiverse" 
scenario of j!4| . In this paper we will be concerned with 
the large scale effects of ultra-light scalar fields on struc- 
ture formation by computing the resultant matter power 
spectrum in a cosmology where a fraction, / = Cl a /Q m , 
of the CDM is made up of a such field. 

If ALPs exist in the high energy completion of the 
standard model of particle physics, and are stable on 
cosmological time scales, then regardless of the specifics 
of the model Tegmark et al have argued [23] that on 
general statistical grounds we indeed expect a scenario 
where they make up an order one fraction of the CDM, 
alongside the standard WIMP candidate of the lightest 
supersymmetric particle. However it must be noted that 
there are objections due to Mack and Steinhardt [25] [26] 
when we consider a population of light fields in the con- 
text of inflation. The problem with these objections is 
that they make some assumptions about what we mean 
by "fine tuning" of fundamental physical theories, which 
is also related to the problem of finding a measure on 
the landscape of string theory and inflation models (see, 
for example, [37]), the so called "Goldilocks Problem". 
Addressing these arguments in any detail is beyond the 
scope of this paper, but we consider the issue sufficiently 
unresolved, and ultra-light scalar fields to be sufficiently 
well motivated as dark matter candidates otherwise, to 
press on regardless in search of phenomenology. 

In the context of generalized dark matter [25] we can 
see the effect of the Compton scale of these fields through 
the fluid dynamics of the classical field. The sound speed 
of a field with momentum k and mass m at a time where 
the scale factor of the FLRW metric is a is given by: 



4m 2 a 2 ' 



k < 2ma 



1; k > 2ma 



(4) 



On large scales the pressure becomes negligible, the 
sound speed goes to zero and the field behaves as ordi- 
nary dust CDM and will collapse under gravity to form 
structure. However on small scales, set by A c , the field 
becomes relativistic and the particles free-stream, sup- 
pressing the formation of structure. This observation is 
our main point of departure to consider the effect of ultra- 
light scalar fields on the matter power spectrum. 

The paper is organised as follows: in Section [TT] we set 
out our system of equations for the cosmological expan- 
sion and perturbations. In Section [nT] we use analytic ap- 
proximations and solutions to identify important scales in 
the evolution of the background field and discuss in more 
detail the issues concerning the fraction of dark matter 
from ALPs, and its production. We go on in Section [TV] 
to set up the initial conditions for perturbations in the 
metric, radiation, dust matter, and the scalar field. We 



then work analytically to approximate the effect of an 
ALP dark matter component on the matter power spec- 
trum and see the emergence of the free-streaming scale 
in the problem. In Section [V] we briefly present numer- 
ical results for the background evolution and the effect 
of ultra-light scalar fields on the cosmic expansion. Our 
main results on the effect of ultra-light scalar fields on the 
matter power spectrum are given in Section[VlJ where we 
give some useful parameterisations to describe them. Fi- 
nally in Section |VII| we discuss these effects and their 
possible detection and outline future directions of work. 



II. THE FORMALISM: EQUATIONS OF 
MOTION 

We work in first order cosmological perturbation 
theory of the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric, in the synchronous gauge, as presented 
in [29|. The line element is: 



ds 2 = a 2 (r)[-rfr 2 + (5 y + h i j)dx i dx :i ] 



(5) 



where a(r) is the scale factor, and r is conformal time. 
The scalar modes of hij can be written as a Fourier inte- 
gral in terms of the two fields h(k,r) and r)(k,r): 



h i:j (x,T) 



d A ke 



i l. 'r 



[kikjh(k, t) + (kikj - -5^)6^,, r)] 

(6) 

where fcj is a unit vector in the direction. 

For a perfect fluid of energy density p and pressure P 
the energy momentum tensor is given by: 

T\ = -(p + 5p) 
T\ = {p + P) Vl 
T j = (P + SP)5) 

(7) 

where p and P are the average density and pressure, and 
5p and 5P represent first order perturbations about ho- 
mogeneity and isotropy. To zeroth order, the Einstein 
equations give the Friedmann equation: 



U 2 = 



8irG 



-a p 



(8) 



where an overdot denotes a derivative with respect to 
conformal time r. The first order equations are: 



k 2 r ] ~ l -Wi= l -a 2 8T\ 
h + 2H- 2k' \ = -a 2 8T l i 



(9) 
(10) 



from which r\ can be eliminated, leaving us with a second 
order equation for h: 



h + Hh = a 2 [6T° Q -5T l 



(11) 
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To couple a scalar field to these equations we compute 
the energy momentum tensor from the potential in the 
usual way: 



V ;Q + 2V)81 



(12) 



Working to first order in perturbations about a homoge- 
neous field: 

0(A j t) = o (t)+&(*,t) (13) 
we have, for a quadratic potential V{4>) = (1 / '2)m 2 <fi 2 : 



Pa =-^-<t>o 
5p a =a~ 2 4>o4>i + ™ 2 ^o0i 



m ^2 



P, =- 



-2 ;,2, 



SP a =a -<p (pi - m-<po<p\ 



{p + P)0 a =aT 2 k 2 ^ 



(14) 

(15) 

(16) 

(17) 
(18) 



The advantage of the synchronous gauge is that all of 
these quantities are independent of the metric perturba- 
tions. 

Next we require the equations of motion for <f>o and 
(f>i, which are found from the Lagrangian, Eqn. [2] (equiv- 
alcntly we could use the conservation equations and Ein- 
stein equations to work directly with the fluid dynamical 
variables of Eqns. 14 15 16, 17] [28 , but making this 
computationally tractable for scalar fields would require 
us to make further approximations that do not always 
hold in the regions of parameter space we are interested 
in): 



O + 2H(f> + m 2 a 2 (/) Q = 
ii + 2H4>! + (m 2 a 2 + k 2 



i = ~^4>oh 



(19) 
(20) 



To obtain the evolution equations for perturbations in 
the dust CDM and the radiation we use conservation of 
energy momentum, T^".^ — 0, in fc-space to obtain the 
first order conservation equations: 



S = -(l + w) + 



h 



= -H{1 - 3w)6- 



w 



SP 



1 + w 1 + w 
The variables 9 and a are defined as: 
(p + P)e = ik j 6T° j 



(21) 
(22) 

(23) 
(24) 



57 j is the traceless component of T l a , a perturbation 
we henceforth ignore; w = P/ p is the equation of state, 
and 5 = Sp/p is the overdensity. 




Working with no baryons coupled to the photon fluid, 
using the CDM particles as the comoving fluid that de- 
fines the synchronous gauge i.e. C = 0, and noting the 
equations of state w y = (P/p)^ = (5P/Sp) 7 — 1/3, and 
w c = (P/ p)c(SP/Sp) c = we have: 



(25) 
(26) 
(27) 



Eqn. [25] can easily be integrated once we have the initial 
conditions to give S c = —\h 1 so that the evolution of the 
matter becomes trivial and we need only work with h. 



III. BACKGROUND EVOLUTION AND 
PRODUCTION OF ALPS IN THE EARLY 
UNIVERSE 

We are interested in scenarios containing a fraction of 
the total energy density today in an ultra-light scalar 
field, therefore we would like to be able to specify f2 a 
in terms of the initial displacement of the field, 4n, or 
equivalently the initial misalignment angle, 0j. To do 
this we look for an analytic solution to the equation of 
motion Eqn. |19| This is most easily done in physical 
time, defined by dt — a(r)dr. In this subsection only, 
overdots will denote derivatives with respect to t, so that 
the Hubble parameter is given by H(t) — a/a. 

We work in reduced Planck units l/m^ = 8irG = 1. 
We rescale to use dimensionless variables t —> Hot, 
H — > H/Hq, (p — > <p/mpi, to — > m/H , where H is Hub- 
ble today, and remain in these variables until we discuss 



the matter power spectrum in Section VI The equations 
governing the background become: 



4> + ?>H<t> a + 777 ^0 = 

tt2 _ Pa(t) n n, 

where the density in ALPs is now given by: 
Pa(t) = + \ m2( t>l 



(28) 



(29) 



(30) 



Eqn. 28 can be solved in terms of Bessel functions if 
we take the ansatz a(t) oc t p 7 which is true in both radia- 
tion dominated (early time) , and matter dominated (late 
time) eras, giving: 

<f>o(t) = a(t)- 3/2 {mt) 1/2 (AJ n {mt) + BY n (mt)) (31) 



with n = ( 1/2) a/ 9p 2 — 6p + 1. We ignore the Y n solution 
since it is singular at early times where we know that 4>q 
should take its value from the misalignment angle. The 
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asymptotic forms of J n tell us how the energy density in a 
scalar field redshifts at early and late times and exhibits a 
well know feature of scalar field evolution in an expanding 
universe. For nit <C 1: 



J n (mt) 



1 



r(n + 1) 



mt 
~2~ 



(32) 



Substituting into Eqn. |31[ along with a oc t p yields: 

cf>o(t)<xt-i p th n (33) 

which gives </>o = const, for both the radiation domi- 
nated era (p = 1/2, n = 1/4) and the matter dominated 
era (p = 2/3, n = 1/2). This in turn shows that the en- 
ergy density remains a constant in this regime: at early 
times the energy density in ALPs redshifts like a cosmo- 
logical constant. 

Later, such that mt 3> 1, we have that: 



J n (mt) 



irmt 



1/2 



cos mt — 



(34) 



Substituting H = p/t now gives: 



fi a (t) 



A 2 m 2 1 
37r a 3 



9 p 2 f n7r 7r 

— t rrr COS mt - 

4 (mi) 2 V 2 4 



1 



oc — +0{(mt)- z ) 



(35) 



for all values of p. We see that at late times the energy 
density in axions redshifts like ordinary matter. 

What these simple observations do not tell us about is 
the transition from cosmological constant (A) behaviour 
to Dark Matter (DM) behaviour, and how this transition 
can affect the expansion rate and age of the universe if 
it contains a significant fraction of DM in ALPs. As 
we will see later there are novel effects even here in the 
background. 

For now we will continue to work analytically and de- 
lineate two important scales in the evolution, an impor- 
tant region of ALP parameter space, and set the initial 
condition on <po for a given f2 a . 

The axion field starts oscillating in the crossover be- 
tween the two asymptotic expressions for the Bessel 
function, when mt osc rs 1, which is the same order as 
the time when the mass overcomes the Hubble drag, 
m m 3H(t osc ). This defines one scale in the problem. 
The background evolution will depend on whether this 
occurs in the radiation or matter dominated era, defined 
by p m (a eq ) = p 7 (a eg ), where p m = p c + p a , the total 
density in matter. If the field has begun oscillations 
in the radiation dominated era then it will be rcdshift- 
ing like matter and contribute as usual when deriving 
a eq ~ f2 7 /fl m . However, if the ALPs begin oscillating in 
the matter dominated era, they will be redshifting as a 
cosmological constant at equality and will contribute neg- 
ligibly to p m . In particular, we can ignore p a (t)a /p c (to), 



so that if the ALPs make up a fraction, /, of the total 
density in matter today, we obtain the modified formula 
for the scale factor at equality: 



1 



On (1 - /) 



(36) 



For ALPs that begin oscillations in the matter era, this 
change to the redshift of equality will have knock-on ef- 
fects for the estimation of other cosmological parameters 
and could possibly place tight constraints on such a light 
species making up an order one fraction of the total dark 
matter density. 

The temperature of the CMB fixes ft 7 ~ 8 x 1CT 5 . 
Then using Vt m ~ Vl c ~ 10 4 £! 7 , simple substitution gives 
that m ~ 10 6 separates fields that begin oscillations dur- 
ing the radiation and matter dominated eras. 

We can estimate the contribution to the critical density 
today coming from ALPs by assuming an instantaneous 
transition from A to DM behaviour and redshifting the 
initial constant energy density from a(t osc ) to a(to) = 1 
as if it were ordinary CDM. We have that: 



1/6 , 1 x 1/2 

t ) \mt Q J 

2/3 

-i.i ; in :. id' 



> io 6 



\mt J 



(37) 



which leads to: 



MU) 2 ; m>10 6 



2 

I) I 



< io 6 



(38) 
(39) 



These expressions can be easily inverted to find an ex- 
pression for the initial condition ^>o(ti,Q, a ). In our code 
we supplement these with an iteratively improved con- 
stant to take into account the ALP effects on t , and 

There are two remaining quantities to be determined, if 
we take the initial scale factor, a, as an input parameter: 
ti and to- We begin in the radiation dominated era, so 
that: 



t, = 



-1/3 



aH 



(40) 



We know to for a matter dominated universe: to = 2/3. 
If we assume for now that ALP effects on the age of the 
universe away from what we expect from a pure matter 
or A dominated universe are small (we will address this in 
more detail in Section|V]) then we can include the effect of 
A on what initial conditions we will need by integrating 



Eqn. 29 The calculation can be found in the standard 
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textbooks [501151]: 
2 



tc 




(41) 



where Q m = Q c + Cl a , which is a monotonically increas- 
ing function of f^A- a cosmological constant makes the 
universe older. 

The difference in our approach here for computing f2 a 
from that in previous works that have been mainly con- 
cerned with the QCD axion, e.g. [11] [32], is that we 
assume no temperature variation to the axion mass, or 
rather we assume that it has reached its zero tempera- 
ture value quickly, and crucially before oscillations of the 
field begin, which since the fields are so light will be at a 
low temperature any how. We consider this a reasonable 
simplification because we do not in general know the tem- 
perature dependence of the mass for a string axion since 
we do not know what instantons will make the dominant 
contribution to the potential. 

If we wanted to restrict our analysis to true axions, 
rather than the more general case of ultra-light scalar 
field ALPs, there is an important region of axion param- 
eter space, known as the "anthropic boundary" , which 
is instructive to locate. True axions are periodic in 
8 = <f>/f a and so the initial misalignment angle has a 
"maximum" at 6^ = tt. In addition, for axions that be- 
gin oscillations in the radiation dominated era, fl a de- 
pends on the axion mass (Eqn. |38[ ). The result is that 
for masses m < 10 12 it is impossible, without taking ac- 
count of anharmonic terms in the potential and tuning 
the initial misalignment arbitrarily close to 7r [11] , for 
axions to produce ft a > 1 and overdose the universe. 
This somewhat alleviates fine tuning problems for these 
light, high-/ a ALPs that lead, in the usual case, to one 
having to tune 0(tj) arbitrarily close to zero to prevent 
overclosure of the universe and other cosmological prob- 
lems [TO] US] . Therefore with ultra-light ALPs the fine 
tuning arguments of Mack and Steinhardt [35] [55] lose 
some of their power. 

The flip side to this is that without tuning the ini- 
tial misalignment arbitrarily close to n and including 
anharmonic effects in the potential it is impossible for 
axions with masses much below the anthropic boundary 
to constitute an order one fraction of the dark matter. 
For example, axions that begin oscillating in the matter 
dominated era, such that tt a is independent of the mass 
(Eqn. [39] )) the maximum possible contribution to the en- 
ergy budget is fl a ~ 4 x 10~ 4 . This observation seems 
to cause problems for the arguments in |24j that any ax- 
ion should contribute an order one fraction of the dark 
matter. 

These conflicting observations on fine tuning for axions 
below the anthropic boundary serve as further motiva- 
tion for the discussion and subsequent decision given in 
Section [TJ to generalise to ultra-light scalar fields with a 
quadratic potential that do not respect the shift symme- 
try of axions when considering such low masses as are 
relevant for the FCDM scenario. 



IV. EVOLUTION OF THE PERTRUBATIONS 
A. Initial Conditions 

To solve for the evolution of the density perturbations 
and compute the resulting matter power spectrum we 
need to find the appropriate initial conditions for the 
perturbations in the various fluid components. 

The initial fluid perturbations can be separated into 
adiabatic and isocurvature components. Whilst axions 
are a source of isocurvature perturbations, and these play 
an important role in constraining axion models by their 
effect on the CMB (see, for example, [HI [55]), we will not 
consider them here, since they do not have a considerable 
bearing on the matter power spectrum. 

Adiabatic perturbations are laid down in a scale in- 
variant way after inflation, and occur in all fluid com- 
ponents from their coupling to the metric perturbation, 
h. Working to lowest order in fcr the coupled equations 
Eqns. [TTJ [26] and [27] can be solved analytically and the 
dominant late time growing mode solution is |29| : 



h = C{k){krf 



-~C(fe)(fcr) 2 
-^C(fc)(fcV) 



(42) 
(43) 

(44) 



where C(k) is fixed by the primordial power spectrum. 
The k dependence be fixed by assuming scale invariance, 
which requires S ~ k 1 ^ 2 and so C(k) = Ck~ 3 / 2 . 

To find the initial condition on cf>i we use the condition 
of zero entropy relating adiabatic perturbations in two 
fluids a and b: 



6 a 



Sab — 1 . 

1 + W a 1 

Sab = S a b = 



(45) 
(46) 



From this one finds that the initial values of </>i and </>i 
both depend on ^o,t(0), the initial value of the derivative 
of 0o with respect to physical time (33] . For ultra-light 
scalar fields that are the PGBs of a spontaneously broken 
symmetry this derivative is zero, since the field is frozen 
at the initial misalignment by Hubble drag, and there is 
no initial velocity. Therefore the initial conditions are 
simply: 



h =0 
h =0 



(47) 
(48) 



The initial conditions derived in |29j are for a numeri- 
cal integration beginning in the radiation dominated era, 
where we also begin our simulations. We can integrate 
the background expansion, Eqn. [8] and find: 



(49) 



G 



which is the final input required to fix the initial condi- 
tions completely, up to the constant C that sets the size 
of the initial perturbations, which we need not specify 
since we will only be considering ratios of power spectra. 



B. Suppression of Structure Formation 

Once the scalar field is deep into its oscillatory phase 
and the background evolution is well described by pure 
matter or radiation domination it is possible to solve the 
equations of motion with a WKB approximation. This 
gives: 



)(t) = A Q 



a(T oa c) 



3/2 



cos 



m / dra(r) 



(50) 



. ,' a{r osc )\ 3/2 
01 (r) = A t [ ] ^<>s 



dr(m 2 a(r) 2 + k 2 ) 1 ' 2 
(51) 

These solutions can then be averaged over their rapid os- 
cillations and an expression for the sound speed [181 128] 
is derived in two asymptotic regimes, as seen earlier in 
Eqn.[4j This momentum dependent sound speed leads to 
the emergence of a new scale in the scalar field evolution: 
kn = ma. For k < kR the sound speed in the density 
perturbations is small and the scalar field behaves as or- 
dinary dust CDM. For k > kR the sound speed in the 
density perturbations goes relativistic and the particles 
free-stream. This similarly defines a time tr after which 
a given mode ceases to behave relativistically. 

A mode of wavenumber k crosses the horizon when 
k ~ Ha, or equivalently kr rs 1. For a given k this 
defines a time of horizon crossing r c . Whether the free- 
streaming scale leads to a suppression of structure, and 
corresponding step-like feature in the matter power spec- 
trum depends on the ordering of the times r c and tr. If a 
given mode enters the horizon once it has already become 
non-relativistic, tr < r c , then the density perturbations 
in that mode will behave as ordinary CDM and there 
will be no suppression of structure relative to the stan- 
dard model. However if a mode is relativistic when it 
enters the horizon, so that r c < tr then structure will 
be suppressed on that scale. We estimate the scale k m 
at which suppression of structure formation begins for a 
field of a given mass, m, at a given redshift, z. 

Smaller and smaller k values arc entering the horizon at 
all times. The mode that entered the horizon at matter- 
radiation equality corresponds to k c {z eq ) s» 0.03ft- Mpc -1 . 
As the scale factor increases so does the boundary value 
for relativistic modes, kR. We require kR < k c , for sup- 
pression of structure at a given redshift z. The large 
scales that we are interested in for the FCDM scenario 
suppress structure formation in modes that entered dur- 
ing matter domination. During matter domination a t 2 , 
which leads to a prediction for the mass dependence of 
the scale k m : 



For heavier scalar fields that suppress structure formation 
in modes that entered during the radiation dominated era 
we have: 



,1/2 



(53) 



so that masses in the range 10 < m < 10 separate 
the regions. However we also know that only masses 
m > 10 6 were oscillating in the radiation era, and there- 
fore it is only for these masses that the derivation of 
kR holds, so that in addition we expect some numeri- 
cal corrections and z dependence to be introduced into 
the expressions above by the transitionary dynamics of 
the background expansion between matter and radiation 
domination, which will effect the expansion rate used to 
derive k m , and due to the background scalar field tran- 
sition between DM and A behaviour. This will be most 
severe for fields that arc still undergoing their transition 
at the redshift of observation, z a i, s . 

How much suppression of structure do we expect rel- 
ative to ordinary CDM? The matter power spectrum is 
given by P(k) = 5^, where 6 m is the total overdensity in 
matter: S m — (5p c + 5p a )/p m - For our purposes it will 
be useful to normalise the power spectrum to one on the 
largest scales. After matter radiation equality, density 
perturbations in ordinary CDM grow like S ~ a. It is a 
well known result [34] that if a fraction f(z) of the matter 
is unable to cluster then perturbations grow as S ~ a q , 
where q= l/4(— 1 + ^25 - 24/(z)). The deviation from 
q = 1 at a given redshift will therefore start at the scale 
k m and saturate at the Jeans scale kj = a\J Hm: this 
is why we expect to see the emergence of "steps" in the 
power spectrum relative to ordinary CDM [TH [35] . 

We can estimate the suppression of power, S, in a step 
using a parameterisation found in [35j found by taking 
the ratio of the two different growth rates of the density 
perturbations: 



2(1-9) 



(54) 



max [a 



osci ^eq 



). The size of a step is then 



k m (m) 



,1/3 



(52) 



where a s 
given by 1 — S. 

An important difference between this work and the 
work in [35j is that we make no approximations for the 
evolution of the scalar field when actually computing 
P(k). 



C. The Scales Involved 

Here we summarise the previous sections by restating 
the important scales to consider when thinking about the 
effects of ultra-light scalar fields on structure formation. 

• A scalar field receives an initial value after symme- 
try breaking and at early times it remains frozen at 
this value by the Hubble drag. A frozen scalar field 
behaves as a cosmological constant; once it begins 
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oscillating it will behave as matter. A field begins 
oscillating when: 

H(t) < m 

• Do oscillations begin in the radiation or matter 
dominated era? We will mostly be interested in 
ultra-light fields that begin oscillations in the mat- 
ter dominated era: 

m < 1CT 27 eV 




— n = o.i 




!!„ = 0.2 




— n« = 0.3 




!1„ = 0.4 





• The energy density today in such an ultra-light field 
depends on its initial value as: 



a, 



• Perturbations in the scalar field have a scale depen- 
dent sound speed, so we can ask: are the perturba- 
tions on a given scale at a given time relativistic? 
The scale k R = ma(t) separates the two regimes. 
On small scales: 

k > k R 

the sound speed is relativistic and the particles free- 
stream, suppressing structure formation. 

• Time dependence of the free-streaming scale and 
the finite size of the horizon mean that suppres- 
sion of structure formation will accumulate on 
scales larger than the free-streaming scale. For the 
ultra-light fields under consideration, suppression 
of structure begins at a scale: 



/ m \ 
V 10- 33 eV/ 



V3 /lOOkms- 



/iMpc" 1 



• The steps in the power spectrum caused by this 
suppression of structure depend on the fraction, /, 
of matter in ultra-light fields. The amount of sup- 
pression can be estimated as: 



S{a) 



&OSC 

a 



2(l-l/4(-l + V25-24/)) 



As one would expect, a larger / gives rise to greater 
suppression of structure, as do lighter fields that 
free-stream on larger scales. 



V. RESULTS IN THE BACKGROUND 

Firstly we show a representative figure, Fig. [T] for the 
evolution of the field (f>o ■ This shows how the initial mis- 
alignment depends on Q a . 

In Fig. [2] we show a representative fit to the numeri- 
cal solutions using the analytic results obtained in Sec- 
tion |III| The fit is made by fixing the field to its initial 



FIG. 1: (j>o versus a for m = 10 3 , various S7 C , S7a = 



0.5 



-0.5 




FIG. 2: A fit for 4> versus a for m = 10 5 , n c = 0.8, fi A = 
0. The solid line is the result of numerically integrating the 
equations of motion, whilst the dotted line is the analytic fit 
of Eqn. 131] 



value before a osc , and then applying the analytic solution 
of Eqn. [3l] after a osc . The analytic solution captures the 
decay envelope and transition well, and the estimate for 
a osc is also a good one. 

There are two important scales in the background evo- 
lution: the redshift of matter radiation equality, z eq , and 
the redshift at which axion oscillations begin, z osc . Both 
of these can be identified on a plot of log a versus log p, 
Fig. [3] and agree well with the expected values, so that for 
example with m = 10 3 the oscillations begin in the mat- 
ter dominated era. Fig. [3] again demonstrates the main 
features of the scalar field evolution, this time through 
the redshifting of the energy density: there is a period of 
constant energy density, and a period where the energy 
density redshifts with the matter. However there is also a 
significant region between a ~ 10 -2 and a ~ 10 _1 where 
the scalar field undergoes a transition in its behaviour. 

One way of looking at the length and significance of 
this transition is to see how it effects the expansion rate. 
In a pure matter, or pure radiation dominated era the 
scale factor evolves as a ~ t p , with a slow, smooth tran- 
sition between the two regimes (see Fig. [4]). Therefore 
plotting H.T versus scale factor will extract the time evo- 
lution of p. When the scalar field begins oscillation there 





logpc 




logp- 
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FIG. 3: logp versus log a for m = 10 3 , fi c = 0.8, fi A = 
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FIG. 4: p versus a for m = 10 , fi c = 0.8, Qa = 0. Insert: 
ALPs + CDM divided by standard CDM alone 



is a more rapid time dependence introduced to p, how- 
ever one can demonstrate that for the ultra-light scalar 
fields of interest here, where oscillations begin around 
t 0(1), that the product prlnr remains small com- 
pared to p and therefore the expansion is still well de- 
scribed by plotting Wt versus scale factor. A plot of the 
exponent p during the transition in the scalar field be- 
haviour estimated in this way is also shown in Fig. [4] 
The relative change in the numerical value of p is small 
during this transition. The transition here begins well 
into the matter epoch, and lasts for approximately one 
order of magnitude in scale factor growth, much shorter 
than the transition from pure radiation to pure matter 
expansion, which lasts over four orders of magnitude for 
the case at hand. Results for different masses of scalar 
field show that the length of the transition as one order 
of magnitude in scale factor is approximately mass inde- 
pendent. These two observations of a relatively short and 
small effect on the expansion rate due to the presence of 
an ultra-light scalar field lead to the prediction that such 
a dark matter component will have a correspondingly 
small effect on the age of the universe, and this intuition 
is indeed borne out in the numerical simulation. 

When a non-zero cosmological constant is included the 
effect on the evolution of the background field is minimal. 



An accelerated expansion rate at late times is included 
simply by altering 4>o(U) through the effect on to as dis- 
cussed in Section |III| a larger t Q will simply increase the 
required initial misalignment for a fixed fraction in ALPs; 
the accelerated expansion dilutes the ALP density. Fur- 
thermore, the mass range of ALPs that will contribute to 
the matter fraction, fi m , today must be deep in their os- 
cillatory regime. Those ALPs light enough that we would 
expect oscillations to begin in a A dominated era cannot 
fulfil this requirement and are overdamped by the in- 
creased expansion rate, as such we will not be interested 
in them since they will not be contributing to the matter 
power spectrum and will simply be adding on to increase 
the effective value of A: the lightest fields we will con- 
sider are of mass m > 10 2 . We conclude that the overall 
effect of A is exactly as for standard ACDM and is unal- 
tered by the presence of an ALP component in the dark 
matter. 



VI. THE MATTER POWER SPECTRUM 

In this section we return to using physical, dimension- 
ful variables. 

In this work, in order to have the smallest and simplest 
set of fits to make to our results, which will be most useful 
for comparison with observations, we will simply fit for 
k m at the centre of any step in the power spectrum, and 
take the saturation at kj as given. We will be interested 
in the function: 



T 2 {k) 



-P^ALPs + CDM 



P ( fc )cDM 
from which we can define the step size: 

S = T (k-max 3* k m ) 



(55) 



(56) 



where k ma x is the smallest scale of interest in the numer- 
ical or observational situation at hand. The dependence 
on m, tt a , and z is implicit. Later we will consider the 
effect of non-zero J7a- 

In Fig. [5] we show T 2 (k) for m = 10~ 29 eV for vari- 
ous values of £l a at z = 0. To fit for S we use the pa- 

for T 2 (k) taken at the largest 



54 



rameterisation of Eqn. 

k = 3 x 10 2 /iMpc -1 in our numerical results, modified 
with the addition of two exponents, fii and fi%: 



S(z) 



{i + zY 



l +(9l v 2(1-8) 



(1 



(57) 



This fit applies when z osc < z eq and when z < zj, where 
zj is the redshift at which kj(z) = 3 x 10 2 /iMpc^ 1 . The 
factor of q contains the dependence on f(z), and z osc 
contains the dependence on the mass. The parameters 
/3i and $2 are chosen by trial and error, and take into 
consideration the approximations in the free-streaming 
scale and oscillation scale used to derive the fit. For those 
fields still far from matter- like behaviour at a given z, we 
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— n„ = o.i 

n a = 0.2 

flo = 0.3 








\ 







Mode k (fcMpc -1 ) 

FIG. 5: T 2 (k) for m = 10" 29 eV, fi A = 




FIG. 6: S versus f(z = 0) = /o for m = 10" 30 eV, 10 _29 eV, 
10 _28 eV, Q A = 0. 



expect a poorer fit. We also expect a poorer fit for those 
fields where z osc ~ z eq . 

In Fig. [6] we show S(f) at z = for three masses 
of scalar field. The fits are used to determine 02, and 
a reasonable match to within a few percent is obtained 
using the value (3 2 — 0.6. Next, Fig.JTJshows S(f(z)) for 
rn = 10 -29 eV with three values of z. The fits are again 
reasonable for a value of f3\ — 0, across two orders of 
magnitude in mass. Particularly, for m — 10~ 29 eV the 
fit is good up to z ~ 100, but breaks down at z ~ 200, 
which is what we expect since at such high redshift the 
field is still very early in its transitionary regime and has 
yet to complete a full oscillation. 

The final quantity we fit for is k m , which we define by: 



T (k m ) — 



(58) 



The fit should be well described by Eqn. 52 and only de- 
pend weakly on z and Q a outside of the transition regime 
in the scalar field behaviour. For sufficiently small z such 
that the field in question has completed its transition we 
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2 = 60 




- 















0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

/(*) 



FIG. 7: S versus f(z) at z = 10, 40, 60 for m = 10 _29 eV, 
n A = 0. 



use the fit: 



^/o ai (l 



)a am l/3 



(59) 



where /o is the fraction in ALPs at z — 0, A is a constant 
of proportionality, and a\, cti are exponents to be fit for. 

Fig. [8] shows the shape of k m as a function of m at 
z = for three values of £l a . At z = all the masses 
of scalar field under consideration have, for the best 
part, undergone their full transition to matter like be- 
haviour, and the shape is well fitted by Eqn. [59] with 
A = 1.25, ol\ = —0.5. With these values essentially nor- 
malising each curve, the shape follows m 1 / 3 well, until 
10- 28 eV. 



52 



and 



m > 10 eV. The discussion surrounding Eqns. 
[53]tells us that these masses suppress structure formation 
in modes that entered when the background expansion 
rate was transitioning from pure radiation to pure mat- 
ter dominated behaviour. In addition, at the time when 
the associated modes entered the horizon the masses in 
question were still undergoing their own transition from 
cosmological constant to dark matter behaviour, which 
will not only effect the background expansion rate but 
also the growth of scalar field overdensities so that nei- 
ther Eqn. [52] nor Eqn. [53] will hold exactly for the mass 
dependence of k m . This departure from the naive picture 
underlines the importance of full numerical simulation to 
determine the effect of ultra-light scalar fields on the mat- 
ter power spectrum. 

Fig.[9]shows the shape of k m as a function of m at three 
different z values. When a scalar field is transitioning in 
behaviour from cosmological constant to CDM, or indeed 
when it is still behaving as pure cosmological constant, 
then the sound speed of Eqn. [4] no longer applies and the 
steps in the power spectrum cannot be simply described 
by Eqn. [52] At high redshift the lower mass scalar fields 
cause the curve to peel away from the m 1 / 3 shape. The fit 



of Eqn. 59 does not attempt to capture the transitionary 
behaviour, and we only note that for fields with a larger 
mass that have completed the transition to CDM there 
is no remaining z dependence in the location of the step, 
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FIG. 8: k m versus m at z = for £l a = 0.1, f2 a = 0.2, 
fi a = 0.3, Q A = 0. 
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FIG. 9: km versus matz = 0, 10, 20 for fl c = 0.8, f2 A = 0. 



i.e. a>2 = 0. 

Given the step size, and the location of the centre of 
the step it is simple to fit the shape of T 2 (k): 



T 2 (k) 



l + S{k/k m )i 

1 + (fc/fcm) 7 



(60) 



The difference plot for this fit against the steps shown 
m Fig. [§ with 7 = 2 is shown in Fig. [lOj and is good 
within around 5% of the total step size, 1 — S. The fit is 
worst around k m , where we expect the most uncertainty 
in our fits, both from the discussion of the fit to Eqn. [59} 
and from our choice to fit for the centre of the step only, 
ignoring the difference between k m and kj. 

When a non-zero cosmological constant is included cor- 
responding to SI a = 0.7, with no change to the fits, they 
remain good at large k and so we conclude that S has 
little or no dependence on Q\. The fit is very poor for 
so we modify Eqn. p9Uo include one more expo- 
nent, as'- 
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FIG. 10: Difference plot for the analytic fit of Eqn. |60| with 
7 = 2 to the numerical results shown in Fig. [5] 
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FIG. 11: Difference plot for the analytic fit of Eqn. |60| with 
7 = 2, using Eqn. |6l] with 03 = 0.75 to fit k m , to numerically 
generated power spectra for m — 10 _29 eV and J1a = 0.7 at 
z = 0. 



With «3 = 0.4 we return to fits good to around 5 — 15% 



^/o ai (l + ^ 2 (l 



n A ) 



a3 m 1/3 



(61) 



of the total step size, as shown in Fig. 11 

We have not considered power spectra for fields with 
m > 10 -27 eV that begin oscillating in the radiation era, 
and suppress structure formation in modes that entered 
during this era, since their behaviour will be much sim- 
pler: scalar fields with these masses will be much better 
fit by the approximations of |35j . and of our fits, in al- 
most all observational situations. They will not effect the 
epoch of equality, and will only effect the expansion rate 
beyond z eqi and, because they will be well transitioned 
to DM behaviour at all redshifts of interest then the scale 
k m will be extremely well fit by 771 1 / 2 . 



VII. DISCUSSION 

In this paper we have explored the cosmological be- 
haviour of ultra-light scalar fields and found that they 
can have a significant effect on the growth of structure 
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at late times. This is not surprising given that we are 
essentially studying the clustering of very light axions, 
i.e. axions whose mass scale is close to the current cos- 
mological horizon. 

Ultra-light axions have, until now, been difficult to mo- 
tivate theoretically. Most of the focus of previous re- 
search has been on QCD or QCD-like axions for which 
there are well defined production mechanisms in the early 
Universe. There have been some studies of what has been 
dubbed Fuzzy Cold Dark Matter [18] , albeit in a slightly 
different mass regime. But with the rise of the string 
axiverse P3]> [35] it makes sense to loosen our usual as- 
sumptions that tie the axion fractional energy, f2 a , to 
the mass scale, m. This means that we have been gen- 
erous with what we deem allowable for the misalignment 
angle- to consider appreciable f2 a s we must allow large 
misalignment angles, pushing the dynamics of the axion 
field towards the anharmonic part of the axion potential. 

We have studied scalar fields whose evolution is such 
that the length scale for oscillations and clustering is 
greater than that set by the horizon size at radiation- 
matter equality. Hence, the transition from "frozen" 
(cosmological constant) to oscillatory behaviour occurs 
after radiation-matter equality. Two interesting new fea- 
tures arise. First of all, the energy density of the scalar 
field behaves as a dark energy at early times and then 
transits to dust like behaviour at late times- this tran- 
sition happens after equality. Second, and as result of 
the previous point, the redshift of equality is shifted by 
a factor (1 — /). If / is substantial, the affect can be ap- 
preciable. This transition between dark energy to dust 
like behaviour will affect the expansion rate and, unlike 
the standard axion picture, it will not be purely dust 



dominated. This effect can be of order 10% for il a ~ 0.3. 

When we look at the evolution of perturbations and 
their effect on clustering we find that, as was expected, 
there is suppression of power on the smaller scales. The 
transition wave number, k m ~ m 1 / 3 , a mass dependence 
which is specific to these ultra-light fields, i.e. with a 
transition after radiation-matter equality. We have re- 
fined the dependence of the scale at which the suppres- 
sion kicks in as well as the amount of the damping on 
small scales, giving us a better handle on how this effect 
depends on parameters such as m and Q a . 

The next obvious step is to see if these various effects 
can be picked out in cosmological observables. These 
ultra-light fields affect the growth rate of perturbations 
and should be detectable by the standard selection of 
methods: the Integrated Sachs- Wolfe effect in the CMB, 
weak-lensing of the large scale structure and galaxy, clus- 
ter and Lyman-a surveys. An intriguing possibility is 
that these effects may contaminate and bias character- 
istic scales in large scale structure such as the Baryon 
Acoustic Oscillation. We will study these different ef- 
fects in a future publication. 
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